Finite Temperature Dynamics of the Spin 1/2 Bond Alternating Heisenberg 

Antiferromagnetic Chain 
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We present results for the dynamic structure factor of the S = 1/2 bond alternating Heisenberg 
chain over a large range of frequencies and temperatures. Data are obtained from a numerical 
evaluation of thermal averages based on the calculation of all eigenvalues and eigenfunctions for 
chains of up to 20 spins. Interpretation is guided by the exact temperature dependence in the 
noninteracting dimer limit which remains qualitatively valid up to an interdimer exchange A ~ 0.5. 
The temperature induced central peak around zero frequency is clearly identified and aspects of the 
crossover to spin diffusion in its variation from low to high temperatures are discussed. The one- 
magnon peak acquires an asymmetric shape with increasing temperature. The two-magnon peak is 
dominated by the S = 1 bound state which remains well defined up to temperatures of the order of 
J. The variation with temperature and wavevector of the integrated intensity for one-magnon and 
two-magnon scattering and of the central peak are discussed. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Gb, 78.70.Nx 



I. INTRODUCTION 

Low-dimensional gapped quantum antiferromagnets 
have received much interest in recent years both ex- 
perimentally and theoretically. They serve as model 
substances allowing to investigate in detail the ef- 
fects of quantum fluctuations and to test theoretical 
models. A particularly simple class of materials in 
this context consists of an assembly of dimers (two 
strongly coupled spins 1/2) which interact sufficiently 
weakly to guarantee that the dimer gap does not close. 
These materials are characterized by a disordered sin- 
glet ground state and a finite spin gap to triplet ex- 
cited states. Materials of this type occur in nature as 
quasi one-dimensional (ID) (Cu(N03)2-2.5H20)i, quasi 
two-dimensional (BaCuSi20e) 2 and three-dimensional 
(KCuCl 3 and TlCuCl 3 )^ compounds. 

Systems of this type are characterized by both quan- 
tum and thermal fluctuations: Thermal fluctuations are 
controlled by the temperature T and the quantum fluc- 
tuations of interest in dimer systems (i.e. those modify- 
ing the basic quantum mechanics of isolated dimers) are 
governed by the interdimer coupling. So far, most in- 
vestigations have concentrated on the quantum aspects, 
excluding thermal fluctuations by working at zero (the- 
oretically) resp. low (experimentally) temperature. The 
theoretical approach is then restricted to the low-lying 
eigenstates of the quantum hamiltonian, which are ac- 
cessible to an approximate analytical treatment or to nu- 
merical calculations. However, at finite temperature of 
the order of or even large compared to the magnetic ex- 
change energy, all levels participate and more elaborate 
approximations resp. numerical approaches are neces- 
sary. 

Recently increasing interest has developped in the fi- 
nite temperature behavior of such systems and exper- 
imental results for Cu(N0 3 ) 2 -2.5H 2 1 ' 5 and TlCuCl^ 
demonstrate that an understanding of the interplay of 



thermal and quantum fluctuations requires additional 
work on the theoretical and numerical side. Our aim 
here is to investigate for the system of weakly coupled 
dimers described above the dynamical structure factor 
over a wide range of frequencies and temperatures. 

To explore from numerical calculations the content of 
a given hamiltonian for finite temperatures is more in- 
volved than for T = since the Lanczos approach which 
allows to treat reasonably large systems (up to 36 spins 
S = 1/2) does not give more than a few low-lying en- 
ergy levels. Instead, the direct approach to dynamic 
properties at finite temperatures requires full exact di- 
agonalization (FED) to obtain energies and wave func- 
tions; this has been done so far for the XXZ S=l/2 chair. 7 
and the frustrated S=l/2 Heisenberg chain* up to 16 
spins. If FED is restricted to energy levels only, thermo- 
dynamic properties such as the specific heat, susceptibil- 
ity and structure factor S(q) are obtained; this has been 
done for a number of dimer type systems such as the 
dimerized and frustrated S = 1/2 chain 9 and the two-leg 
ladder—. These numerical approaches to finite tempera- 
tures have been supplemented by approximate analytical 
approaches such as including thermal occupation factors 
in the mean field approach&Al and a strong coupling field 
theoretic apporach— 

We describe the one dimensional (ID) 5 = 1/2 chain 
with alternating isotropic antiferromagnetic nearest- 
neigbour exchange (the bond alternating Heisenberg 
chain, BAHC) as a system of N/2 unit cells with two 
spins each and two exchange constants, the intradimer 
exchange J and the interdimer exchange AJ, using the 
following hamiltonian: 

N/2 

), (i) 

We assume J > and apply periodic boundary condi- 
tions. For A = the ground state of the system consists 
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of singlets on the intradimer bonds (n, 1) — (n, 2). These 
local singlets can be excited to triplets which develop 
into a band of gapped excitations when switching on A. 
At higher energies multiparticle excitations dominate the 
spectrum. In the limit A = 1 we arrive at the well known 
Heisenberg antiferromagnetic chain (HAF) with pairs of 
S = 1/2 spinons as lowest gapless excitations. Other 
related models are described in Ref. Il3l 

The dynamics of such a system is most appropriately 
discussed in terms of the dynamic structure factor, i.e. 
the Fourier transform of the time dependent spin cor- 
relation function which, apart from well known pref- 
actors gives the spectral weight measured as the mag- 
netic inelastic neutron scattering (INS) cross sectionM. 
The dynamic structure factor of the BAHC at low tem- 
peratures (T <C J) is characterized by a peak due to 
singlet-triplet transitions as most prominent feature and 
continua of multiparticle excitations at higher energies. 
From extended work on the T — quantum mechanics of 
this system (Refii&aiSiiLiSiiii and references therein) it is 
known that finite interdimer coupling A leads to consid- 
erable modifications in energies and transition strengths 
for one and two magnon processes as well as to the 
emergence of bound states below (and above) the two 
magnon band. Here we present results from FED for up 
to 20 spins for the longitudinal dynamic structure fac- 
tor (DSF) S zz (q, oj) for a wide range of temperatures, 
varying from T < J to T > J. Our discussion will con- 
centrate on the variations of the one and two magnon 
lineshapes with temperature (with particular emphasis 
on the bound state) and on the temperature dependence 
of the contributions from the temperature induced 'cen- 
tral peak' at u> ~ which is due to intraband transitions. 
This peak is the signature of spin diffusion in the classical 
limit; its temperature dependence has not been discussed 
so far in a microscopic context, but it is clearly seen in 
recent INS experiment A 

For N spins localized on sites jq on a ID lattice the 
dynamical structure factor is defined as 

(S?(t)S?(t = 0)) (2) 

The superscripts a, (3 denote the spin components and 
the brackets (• • •) thermal expectation values (which 
for T — reduce to groundstate expectation values 
(0| • • • |0)). For the isotropic hamiltonian (no magnetic 
field) of Eq. 0) the dynamic structure factor is diagonal 
with three equal components. It is therefore sufficient 
to restrict the discussion in the following to one compo- 
nent, <S zz (q, ui). In the presence of a magnetic field, a 
case to be treated in subsequent work, we will have to 
distinguish between S zz (q, oS) and S A (q,w). 

The remainder of this article is organized as follows: In 
section[n]we review and discuss simple limiting situations 
and our numerical approach, in section IIIII we present 
results for the BAHC with A up to 0.3 and temperatures 



up to 4 J and in the final section Hvl we discuss the range 
of validity of our approach and give our conclusions. 



II. LIMITING CASES AND NUMERICAL 
APPROACH 

The BAHC has lattice sites x, = nb + pd/2, n = 
1...N/2, p = ±1. Here d is the intradimer distance, the 
direction of d in general differs from the chain direction 
given by b. At finite temperatures the general expression 
for S zz (q,ui) is written as 



i 

E K^l S k |^>| a *(w-(w fc -Wi)) (3) 

k 

where Z is the partition function and we have introduced 
the Fourier transformed spin operators 

N/2-1 

S z = — Y y e lq (™ b +P d / 2 ) S z . (4) 

q y/N ^ ^ ' P 
v n— p— ±1 

We begin by reviewing results for the two limiting cases 
T = (pure quantum effects) and A = (pure thermal 
fluctuations). The BAHC at T = with its singlet-triplet 
gap is one of the simplest systems to study strong quan- 
tum fluctuations in numerical and analytical approaches. 
This has been done by many different methods, among 
them the random phase approximation^^ and series ex- 
pansions in the coupling strength .^iLiLiSiiS. The appli- 
cation of these methods has concentrated on two types of 
low-lying excitations, i.e. eigenstates of the hamiltonian, 
above the singlet ground state: (i) The basic quantum 
excitation is the first excited triplet, i.e. the propagat- 
ing one magnon excitation (one excited dimer) (ii) two 
magnon excitations which form a continuum with the 
possibility of bound states. These states are decorated 
by admixtures of multimagnon excitations which become 
more and more important when the interdimer coupling 
A increases. These methods have established the exis- 
tence of bound states of two excited dimers below the 
two magnon band, two singlets and two triplets^, for 
wave vectors close to 7r and energies and transition ma- 
trix elements for the one and two magnon bands have 
been obtained up to 13 th order— from series expansions. 
We remark that the analogous calculation including 3D 
interactions has only been done up to 3 rd order 19 . These 
results are in excellent agreement with experiments at 
low temperatures. 
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In the limit of vanishing coupling, A = 0, on the other 
hand, the BAHC reduces to an assembly of noninteract- 
ing dimers and its exact dynamical structure factor can 
be easily calculated for all temperatures. In spite of the 
simplicity of the calculation the result is instructive as a 
guide to the general case. The full result for S zz (q,uj) in 
this limit is 



1 i 



{2(l + cosq-d)e" /3J (5(w) 



4 1 + 3e-/" 

+ (1 - cos q • d) (5(w — J) + e _/3J (5(w + J)) } (5) 

This expression is valid for any number of spins N. 
The wavevector q is in chain direction and q takes 
only the N/2 discrete values q = (2irm)/{Nb/2), m = 
1,2,... N/2 of the lattice of dimers (lattice constant b). 
Eq. (jSJ shows the two types of contributions which sur- 
vive the limit A — > 0, i.e. which are present in the non- 
interacting dimer limit: (i) one magnon excitations with 
energy co = ± J (related in strength by the detailed bal- 
ance factor e _,3J ); (ii) a contribution at to = which 
results from transitions (diagonal terms for S zz ) within 
the excited dimer triplet. It therefore carries a factor 
e~P J and will develop into a 'central peak' for A ^ 0. 
Since all transition processes in this limit are local- 
ized, the only nontrivial wavevector dependence results 
from form factors related to the internal dimer structure: 
the intradimer distance d determines the 'dimer struc- 
ture factor' (1 — cosq • d) and the corresponding factor 
(1 + cosq ■ d) for the central contribution. 

The result of Eq. (J5J is independent of N since for 
independent dimers, A = 0, the correlation length van- 
ishes and the exact result is obtained already for N = 1. 
Increasing N only changes the number of wave vectors 
q and their positions. This guarantees the extensivity 
of the sum rule £ q / du>S zz (q,,u) = N/4 20 . The abso- 
lute intensity of the one 'magnon' peak at e.g. q = 2n 
develops from 1/2 at T = to 1/8 at T -> oo. We 
note that the various contributions to this intensity from 
states with a different number of excited dimers depend 
on N as can be identified easily from Eq. ©: The pref- 
actor in Eq. (JSJ, 



1 _ (l + 3e-^ J ) N ^ 
l + 3e-/"~ (l + Ze-PJ)"/ 2 



2-1 



(6) 



originates from the partition function in the denominator 
and the phase space factor resulting from excited dimers 
which do not participate in the transition in the numera- 
tor. The expansion of the numerator thus gives as a series 
the relative strength of the various contributions with 
energy J starting from a state with M excited dimers: 
The contribution of states with M excited dimers is pro- 
portional to (^ 2 ) (3 e - ^ J ) M . For T — » oo this reduces 
to { N ' 2 ) 3 M , the de generacy of states with M excited 
dimers. From this follows that the various contributions 
to the dynamic structure factor from states with M ex- 
cited dimers have a well defined maximum (for N 3> 1) 




FIG. 1: Density of states for N = 20, A = 0.3, q = -it. Each 
point in the figure represents one multiplet. 



at M = M = zN with z = 3/(e /3J + 3). This means 
M w 0.5247V for T = J and M = 0.757V for T -> oo. 

It should be noted that even for T — > oo the result 
of Eq. lO is a quantum mechanical result reflecting the 
discrete energy spectrum of a S — 1/2 dimer. It does 
not agree with the classical limit which is defined by the 
limit h — ► 0, i.e. for spins which are classical vectors 
or alternatively in the limit S — > oo. From the latter 
approach follows that the frequency regime < ui < 2JS 
fills continuously with excitation strength in the classical 
limit. The result of Eq. (JSJ can clearly also be generalized 
to include an external magnetic field, a case which will 
be treated in a forthcoming paper. 

Allowing for finite coupling A several features appear: 

(i) the contributions at to = and u> = ±J acquire a 
finite width owing to the dispersion of the excitations, 

(ii) additional contributions at lu w ±pj appear in in- 
creasing order in A, (iii) the intradimer form factors get 
modified 19 . For small coupling A <C 1, two essential 
characteristics of the noninteracting limit survive qual- 
itatively as shown for N = 20, A = 0.3 in Fig. [2 The 
correlation length remains small and the density of states 
continues to show a well defined sequence of peaks with 
maxima at E « pj up to the upper bound of the spec- 
trum and with widths of the order of A independent of 
energy. Thus the largest contributions to the transitions 
uj w J (lu w 2 J), which start at energies E w Mo J will 
be included for N > 16 (JV > 20) even for T -> oo. We 
conclude that these system sizes will be sufficient to ob- 
tain reasonable results for the BAHC structure factor at 
moderate A and all temperatures in its most interesting 
frequency regimes. 

For our numerical calculations we have used the House- 
holder algorithm to obtain all eigenvalues and eigenvec- 
tors of the hamiltonian in the subspace of constant S z ot 
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and wavevector. Calculations were performed on work- 
stations with Pentium IV processors at the ITP, Han- 
nover University and on the JUMP supercomputer at 
NIC Jiilich. Calculations for N — 16 are comfortably to 
perform (matrix dimension is about 2000, required mem- 
ory about 100MB). A few calculations have been per- 
formed for N — 20 for one value of A and the wavevec- 
tor and frequency regime where the bound triplet two 
magnon state is found; here memory (a few GB) and re- 
quired CPU time are at the limit what can be reasonably 
managed at present. 



S**{q,u) for N = 16, A = 0.3, T = 0.5 




III. RESULTS AND DISCUSSION 



S zz {q,uj) for N = 16, A = 0.3, T = 2.0 



In this section we present our results for 5 zz (q, u>). 
We cover the full range of temperatures (in units of J), 
T = 0.1.... T = 4 (the latter being representative for 
T — > oo) and wave vectors q (in units of b~ r ) from to 2ir. 
As standard values we use N = 16 and A = 0.3. We find 
that the variation with A for smaller values can be safely 
found from interpolating the results presented here. The 
number of spins strongly influences the level density and 
the discreteness of reciprocal space; however, the posi- 
tions of levels which occur at identical wavevectors for 
different N (e.g. the one magnon state and the 5=1 
bound state below the two magnon band at wavevector 
7r) do not change with TV going from N = 12 to N = 24 
(values for N = 24 were obtained from Lanczos diago- 
nalization). In this work we consider the simplest case 
d = b/2. Then d is the distance between all spins and it 
is sufficient to specify q — |q| in the following, d = b/2 
implies that the effects from the dimer form factors are 
particularly simple: in lowest orders the central contri- 
bution vanishes for q = 2n and the one magnon con- 
tribution for q = 0. Temperatures and frequencies are 
given in units of J throughout this section. The complete 
variation of S zz {q,oj) with wavevector and frequency for 
N = 16 and A = 0.3 is shown in Fig. [21 for two tempera- 
tures, T = 0.5 and T = 2. The small contribution of the 
two- magnon excitations has been enlarged by a factor 20. 
The remaining figures display the details of the influence 
of temperature on the various interesting aspects of the 
spectra. The scale of S zz (q,uj) is set by the value for 
the one-magnon peak for noninteracting dimers at T = 
and q = 2tt, i.e. 1/2 in the units of Eq. J2J). For all 
presentations the transitions strengths have been added 
up within frequency bins of mostly Aw = 0.01 (0.002 in 
Fig. EJ). 

Fig-EJiUustrates the evolution with temperature of the 
one-magnon peak for N = 16 and A = 0.3. With in- 
crease in temperature, the peaks extend over a nearly 
constant range in frequency, although the decrease of 
the maximum of intensity formally implies an increase 
in the width at half maximum. For both q — 7r and 
2ir the position of the maximum shifts to lower frequen- 
cies, but the peak develops a marked asymmetry with 
more intensity on the high frequency side. This effect 




FIG. 2: Overall picture of the dynamic structure factor 
S zz (q,uj) for N = 16, A = 0.3, (a) T = 0.5, (b) T = 2.0. 
The intensity of the two magnon peaks is enlarged by the 
factor 20. 



S"(q, a) tar <V-16, A-0.30, q-n and q =Stt 




r=o.50 
r=i.oo 

T-2-00 



0.6 0.8 1 1.2 1.4 



0.6 0.8 1 1.2 1.4 



FIG. 3: (Color online) One-magnon peak for JV = 16, A = 0.3: 
(a) q = 7r and (b) q = 2n. Temperatures are T = 0.5 (dotted 
lines), T = 1.0 (dash-dotted lines), T = 2.0 (full lines). The 
arrows indicate the frequencies of the one magnon transition 
at T = 0. 
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is most pronounced for q = 2ir, the wave vector of the 
energy gap (and may not be the case at all for small 
values of q, see Fig. |2J). Thus the detailed description 
of the lineshape provides an understanding for the puz- 
zling observation that the gap energy seems to increases 
with temperature^ (comparable to an analogous obser- 
vation in the Haldane chain 21 ). It would be interesting 
to compare the lineshape to the result of the theoretical 
approach of Ref. ^2- However, with N = 16, a contin- 
uous lineshape results only for temperatures above the 
energy gap and is thus in the nonuniversal approach of 
the approach of Ref. 0. When we look at the micro- 
scopic origin of the peak broadening, we find that e.g. at 
T = 1.0 the basic transition from the ground state to the 
one magnon excitation at q = ir is responsible for 30% of 
the weight at that frequency and that for the neighbor- 
ing frequencies transitions starting from the one magnon 
band contribute about 10% of the total weight, whereas 
the by far largest part of the intensity originates from 
transitions starting at states with two or more excited 
dimers. 

Fig. ^illustrates for N — 16, q — 7r/4 and tt the evolu- 
tion of the central peak with temperature (for A = 0.3). 
At low T and A <C 1 the shape of the central peak is 
dominated by transitions inside the weakly interacting 
one magnon band with q— independent matrix elements. 
The physical process may be thought of as an external 
probe accelerating a thermally populated excitation and 
is thus similar to well known processes in soliton bearing 
ID magnets^. As in these models, the limiting form of 
the structure factor is 



N=16, A = 0.30, q=TT/'4 and q=n 



v /l-w 2 /w ra (?) : 



0(o&(g)-a; 2 ) (7) 



with u) m {q) = Asin(g/2) just from phase space effects 
(a small variation with temperature resulting from the 
dispersion has been neglected). Numerically we find 
that the cutoff of the central peak remains localized at 
lu rs Lu m (q) for all wavevectors and temperatures (apart 
from corrections of C(A 2 )). The lineshape, however, can- 
not be expected to be reproduced since only a few dis- 
crete transitions in the one magnon band are available for 
N = 16 and dominate the spectrum at the lowest tem- 
peratures. Nevertheless, for q — tt,T — 0.5 the inverted 
lineshape of Eq. |J7J| starts to become visible even with 
this restriction. Between medium (T = 0.5) and high 
(T = 4) temperatures a crossover of the lineshape from 
square like to Gaussian is observed for small wavevectors. 
In view of the experimental 2 ^ and theoretical 2 ^ 2 ^ discus- 
sion of spin diffusion in gapped ID magnets we have also 
made a rough analysis of the linewidth variation with 
wavevector. We find that the T = 4 spectra fit surpris- 
ingly well to a Gaussian (but not to a Lorentzian) and the 
linewidth increases by a factor of about 5 upon doubling 
the wavevector from q — it/ 4 to n/2. This indicates that 
the crossover to spin diffusion like behavior governs the 
behavior of the DSF in the regime considered here. The 
decrease of the overall intensity with temperature due to 



— 7=0.50 

- - T = 2.00 
T = 4.00 





-0.4-0.2 0.2 0.4 



-0.4-0.2 0.2 0.4 



FIG. 4: (Color online) Central peak for N = 16, A = 0.3: 
variation with temperature up to T — 4 for (a) q — tt/4 and 
(b) q = ty 



the thermal occupation factor will be discussed below, 
see Fig. 

Fig. shows the two-magnon peak at q — tt (where 
the bound states are most clearly visible) and q — 2tt for 
N = 16, A = 0.3 in the temperature range T = 0.5 ... 2 
(note the enhancement in frequency and intensity scales) . 
For q = tt the low temperature spectra are entirely dom- 
inated by the S = 1 bound state at u « 1.938..., whereas 
the continuum (which at T — is smaller by a factor of 
A 2 in intensity^) plays no significant role. The bound 
state remains clearly visible up to T = 1 and then dis- 
appears in parallel with rapid decrease in the integrated 
intensity of the two-magnon peak for temperatures above 
T = 1 (see also Fig. below). This reflects the fact that 
dimers become independent of each other with increas- 
ing temperature such that the correlations between spins 
in two different dimers required for a finite two-magnon 
peak disappear. In Fig. we present a comparison be- 
tween results for N — 16 and N = 20 (available only for 
q = 7r and the limited frequency range 1.84 < w < 2.02). 
In going from N = 16 to N — 20 some improvement 
is obtained, the continuum becomes smoother and the 
bound state less dominant, but the main characteristics 
are unchanged. Thus the increase in N is not really cru- 
cial. The spectra invoke the impression that at interme- 
diate temperatures additional transitions (in particular 
the S = state at w « 1.868 which could be reached 
from a thermally excited S = 1 state) become visible. 
This, however, is misleading as the comparison between 
the data for N = 16 and N = 20 shows: A change in 
the number of spins leads to a different set of allowed 
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N-16, A-0.30, q-n and q-Zjr 




1.7 1.8 1.9 2 2.1 2.2 1.7 1.8 1.9 2 2.1 2.2 



N = 16 and fV=20, \=0.30, q=7T 



5xl0- 5 



T=0.5 



1.85 1.9 1.95 2 



1.85 1.9 1.95 2 



FIG. 5: (Color online) Two magnon peak for N = 16, A = 0.3 
and different temperatures (symbols as in Fig. |5J|: (a) q — n 
(the S — 1, T = 0.1 bound state contribution is at 0.00318 
outside the frame), (b) q = 2-n (the two largest T = 0.1 
contributions are outside the frame at 0.0109 for uj ~ 1.73 
and at 0.00961 for uj « 1.99) 

wave vectors and thus, from the accompanying change in 
the initial and final energies, to trivially different exci- 
tation frequencies (although the energies at e.g. q = tt 
are remarkably independent of N). For q — 2n there is 
no bound state; the two discrete transitions at T ~ 
in Fig. 03 are the remainder of the continuum for the 
limited number of spins 16. 

For a more global analysis we show in Fig. [7| the in- 
tensities I n (q = it, T) obtained from integrating S zz (q — 
tt,lu) over frequency intervals (n — ^) < ui < (n + i). 
The results for q = 7T are representative also for the 
other q~ values. These 'exclusive structure factors' show 
most clearly the relevance of quantum fluctuations with 
increasing interdimer coupling A at low temperatures. 
Quantum fluctuations are most important for the one- 
magnon peak I\ at low temperatures. Ii(q = ir, T) comes 
close to its classical limit independent of A at T w 2, 
whereas the temperature dependence of Io(q = 7r, T) fol- 
lows closely the A = result of Eq. JSJ for all tempera- 
tures. The temperature variation of / 2 (<7 = tt, T 1 ) shows 
clearly that an observation of the two-magnon peak re- 
quires low temperatures T < 1. Transitions correspond- 
ing to exciting more than two dimers are extremely small 
in magnitude: A typical number is Is(q — 2n, T — 0.5) ~ 
2 • 10~ 5 . From further calculations of integrated inten- 
sities we conclude that the number of spins is not very 
important for the overall features: N — 12 gives nearly 
undistinguishable results. 

Fig. [H] shows the influence of temperature and inter- 



FIG. 6: (Color online) Two magnon peak for q — ir,\ — 0.3 
and N — 16 and iV = 20 (for the frequency interval 1.84 < 
uj < 2.02) in comparison: (a) T = 0.5 and (b) T = 1. Here, 
narrow frequency bins, Aoj = 0.002, are used. 
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FIG. 7: Temperature dependence of the 'exclusive structure 
factors' I n (q, T) for N = 16, A = 0.3 and q = n. n = 0,l (left 
scale, lines give the result for A = from Eq. (y) and n = 2 
(right scale) 



dimer interaction on the dimer structure factor (1 — 
cos qd) for the one magnon intensity and the correspond- 
ing factor (1 + cos qd) for the central peak. In the anal- 
ysis of experiments this factor is usually considered as 
a prefactor independent of A and uj, although it is clear 
already from the zero temperature analysis^ that cor- 
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f„(q, T) for N=16, A = 0.3D 




FIG. 8: Wave vector dependence of the 'exclusive structure 
factors' I n (q,T) for (a) T = 0.5 adn T = 2.0. Lines for 
n — 0, 1 give the result for A = from Eq. JSJ. Units for I2 
(right scale) are enhanced by a factor of 20. 

rections are present in higher orders in A. Fig. [S] shows 
an overall agreement between the numerical results for 
A = 0.3 and those of the noninteracting limit although 
closer inspection reveals significant differences: The ra- 
tio h(q = 2n,T = 0.5) /h{q = n,T = 0.5) e.g. is about 
2.5 for A = 0.3 compared to 2 for A = 0. Fig. [S] also 
displays the wavevector dependence of the two magnon 
peak which is absent in the noninteracting limit: The in- 
tensity of this peak depends strongly on wavevector and 
wave vectors q > tt are most favorable for an observation 
of these processes. 

IV. CONCLUSIONS 

We have calculated the dynamic structure factor of the 
bond alternating Heisenberg chain based on full exact di- 
agonalization for chains with up to 20 spins 1/2. This 
allows us to obtain results at finite temperatures which 
cover the complete temperature range from T <C J to 
T» J. Our results are for dimer like chains with suffi- 
ciently small interdimer exchange AJ. We find that the 
characteristics of the noninteracting dimer limit (where 
the full temperature dependent structure factor, Eq. J£J, 
is easily available analytically) describe much of the in- 
teracting BAHC. This is demonstrated here for A = .3 
(but is essentially true up to typically A w 0.5) and ap- 
plies in particular to the existence of the central (zero 



frequency) peak and to the temperature and wave vec- 
tor dependence of the integrated intensities ('exclusive 
structure factors') for the central and the one magnon 
peak. This large range of validity of the dimer picture 
is in agreement with the observation^ that the dimer 
picture gives a good account of excitation strengths up 
to and including the isotropic limit of the HAF at zero 
temperature. 

Experimental and theoretical interest in the dynamics 
of the BAHC results from the possibility that this theo- 
retical model and its realizations in a number of materials 
might serve as simple model systems to discuss the in- 
terplay between temperature and quantum fluctuations. 
Our results reveal that this interplay becomes apparent 
in a number of points which are accessible to experimen- 
tal, in particular neutron scattering, observations: 

(i) The shape of the central peak which is a prominent 
feature of the BAHC, displays the crossover with tem- 
perature from the noninteracting particle like behavior at 
low temperatures to the diffusive behavior at high tem- 
peratures. 

(ii) The one magnon peak develops an asymmetric 
linewidth (with appreciable strength on the high fre- 
quency side) with temperature. This is particularly evi- 
dent for wavevector q = 2n at temperatures T > 1 and 
appears to describe in more detail the upward shift in 
gap energy with temperature noted in approximate the- 
oretical approaches. 

(iii) The two magnon peak (around q — n) is dominated 
by the bound triplet state on top of a small continuum 
which is smooth at all temperatures and wavevectors. 

The results as presented here are restricted to the 
basic situation of equidistant spins and isotropy implying 
the absence of an external magnetic field. In subsequent 
work we will extend the calculations to cover the finite 
temperature dynamics in an external magnetic field and 
also to describe real materials such as Cu(NOs)2 and 
(VO)2P20y. We will also investigate in more detail the 
regime A > 0.5, approaching the isotropic (Heisenberg) 
limit in order to study stronger deviations from the 
dimer like situation . 
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